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Abstract 

The kinematic wave model of traffic flow on a road network is a system of hyperbolic 
conservation laws, for which the Riemann solver is of physical, analytical, and numerical 
importance. In this paper, we present a Riemann solver at a general network junction. In the 
Riemann solver, we replace the entropy condition in [25] by a local, discrete flux function used 
in Cell Transmission Model [11]. To enable such an entropy condition, which is consistent 
with fair merging and first-in-first-out diverging rules, we enlarge the weak solution space by 
introducing interior states on a set of measure zero, associated with stationary discontinuities 
at the junction. In the demand-supply space, we demonstrate that the Riemann problem is 
uniquely solved, in the sense that stationary states and, therefore, kinematic waves on all links 
can be uniquely determined from feasible conditions on both stationary and interior states as 
well as the entropy condition that prescribes boundary fluxes from interior states. In addition, 
the resulting global flux function is the same as the local one. Thus the flux function is both 
invariant and Godunov. 

Key words: Lighthill-Whitham-Richards model; kinematic waves; Riemann problem; supply- 
demand space; stationary states; interior states; entropy conditions; fair merging; First-In-First-Out 
diverging; demand level; supply level; critical demand level; Cell Transmission Model 

1 Introduction 

A better understanding of traffic dynamics on a road network is critical for improving safety, 
mobility, and environmental impacts of modern surface transportation systems [45]: practically, 
it is helpful for efficient implementations of ramp metering [41], evacuation [46], signal control, 
and other management and control strategies; theoretically, it can yield better network loading 
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Figure 1 : A grid network 



models for many other studies [49]. In a road network, e.g., a grid network shown in Figure 1, 
vehicular traffic dynamics can be described by cellular automata models [39; 10] and car-following 
models [18] of individual vehicles' movements, fluid dynamic models of continuous car-following 
behaviors [42; 48; 2; 50], the Lighthill-Whitham-Richards (LWR) kinematic wave model [37; 43], 
or regional continuum models [4; 24]. The traditional LWR model describes traffic dynamics 
of homogeneous vehicles on a virtually one-lane road as combinations of shock and rarefaction 
waves and can be analyzed with theories of hyperbolic conservation laws [33]. With the right 
balance between physical reality and mathematical tractability, kinematic wave models have been 
successfully extended to study more complicated traffic dynamics of heterogeneous vehicles [5], on 
multi-lane roads [13], and through network junctions [25]. 

In a road network, such bottlenecks as merging, diverging, and general junctions play a critical 
role in initiating, propagating, and dissipating traffic congestion. Some interesting traffic dynamics 
can be caused by interactions among these network bottlenecks: for examples, a beltway network 
can be totally gridlocked [12], and periodic oscillations can occur in a diverge-merge network 
[26]. Thus efforts are warranted to develop both physically realistic and mathematically tractable 
kinematic wave models of network traffic dynamics. Since traffic dynamics inside a link can be 
described by the LWR model and are well understood, the most important component of network 
kinematic wave models is related to how merging and diverging behaviors would impact the 
formation of shock and rarefaction waves at a general network junction shown in Figure 2, which 
has m upstream links and n downstream links. In the literature, there have been three lines of 
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Figure 2: A general network junction 



research into traffic dynamics through a general network junction. In the line of discrete Cell 
Transmission Model (CTM) [11; 34], boundary fluxes through a junction during a time interval are 
prescribed from adjacent cells' conditions based on macroscopic merging and diverging rules. In 
the line of continuous models [25], shock and rarefaction waves on all links are analytically solved 
with jump initial conditions by decoupling the Riemann problem at the junction into m + n Riemann 
problems on individual links. In the third line of continuous models [29], the discrete flux functions 
in CTM are used as decoupling conditions, and it was shown that the Riemann problem can be 
uniquely solved. The first two approaches bear their respective limitations: the CTM approach has 
been verified by empirical observations but cannot be applied to obtain such analytical insights as 
shock and rarefaction waves at a junction; the continuous model by [25] is mathematically tractable, 
but the decoupling method based on an optimization problem is not directly associated with any 
physical merging or diverging rules. In contrast, the third approach integrates the physical merging 
and diverging rules in CTM into the analytical framework of [25] is very promising for further 
studying network traffic dynamics. 

Since the kinematic wave model of network traffic flow is a system of hyperbolic conservation 
laws on a network structure, solutions to the Riemann problem at a junction, in which all links carry 
constant initial conditions, but discontinuities can occur at the junction, are of physical, analytical, 
and numerical importance: physically, they can define physical merging, diverging, and other 
behavioral rules; analytically, a system of hyperbolic conservation laws is well-defined if and only 
if the Riemann problem is uniquely solved [6]; and numerically, they can be incorporated into 
the Godunov finite difference equations [19]. In this study, we present a Riemann solver of the 
network kinematic wave model using the solution framework of [29]. Note that the Riemann solver 
is analytical in the sense of [15], different from numerical ones as discussed in [44]. The new solver 
is based on that in [25]: in the Riemann solutions, a stationary state arises on a link along with a 
shock or rarefaction wave, which is determined by the Riemann problem of the LWR model on the 
link with the initial and stationary states; the stationary state should be inside a feasible domain, 
such that the shock or rarefaction wave propagates backward on an upstream link and forward on 
a downstream link; and the constant in- or out-flux of a link equals the stationary flow-rate. The 
remaining piece in the Riemann solver is to introduce an entropy condition in the sense of [25] 
such that "this condition gives a unique solution at least for Riemann initial data". But different 
from that in [25], the Riemann solver in this study uses a discrete flux function, which is defined in 
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terms of upstream demands and downstream supplies, as an entropy condition. Such a discrete flux 
function is originally developed within the framework of CTM, models conflicts among merging 
and diverging traffic streams at the aggregate level, and is therefore a natural choice as the entropy 
condition to pick out unique physical solutions. To incorporate the new entropy condition, we 
enlarge the function space for weak solutions to the Riemann problem by introducing on each link 
an interior state, which is local and takes no space (of measure zero) right next to the junction. That 
is, the Riemann solution may be neither left- or right-continuous at jump discontinuities, and such a 
function space is still the same as that for traditional weak solutions. Then the entropy condition is 
introduced so that boundary fluxes through the junction be determined locally from interior states 
with numerical CTM flux functions. In addition, based on the equivalence between traffic density 
and the demand-supply pair, we solve the Riemann problem in the demand-supply space and show 
that stationary states exist and are unique. In this study, the local flux function is the global flux 
function derived in [29] . 

The rest of the paper is organized as follows. In Section 2, we present the network kinematic 
wave model and review Holden and Risebro's Riemann solver and discrete CTM flux functions. 
In Section 3, we present a new Riemann solver. In Section 4, we solve the Riemann problem. In 
Section 5, we discuss some further properties of the Riemann solver. In Section 6, we present some 
concluding remarks. 

2 A system of hyperbolic conservation laws at a network junc- 
tion and the Riemann problem 

In this study, we consider traffic dynamics in the junction network shown in Figure 2. We denote the 
set of upstream links by A = { 1 , • ■ • , m} and the set of downstream links by B = {m + 1 , • ■ • , m + n}. 
For link a G A, we introduce link coordinates (a,x a ), where x a G (— °°,0); for link b G B, we 
introduce link coordinates (b,Xb), where x\, G (0,°°). At a point (a,x a ) (a GAUB) and time t, 
we denote the total density, speed, and flow-rate by k a (x a ,t), v a (x a ,t), and q a (x a ,t), respectively. 
Hereafter we omit (x a ,t) from these variables unless necessary. Here we assume that vehicles 
of different paths, classes, or other attributes have the same characteristics, and that each link is 
homogeneous with a location-independent number of lanes, free-flow speeds, curvatures, slopes, 
and so on. Then we have the following fundamental diagram of flow-density and speed-density 
relations [20]: 

V a = V a (k a ), 

q a = Q a (k a ) =k a V a (k a ). 

Generally, Q a (k a ) is a unimodal function in k a and reaches its capacity, C a , when traffic density 
equals the critical density k ajC - If traffic density k a is strictly smaller than, equal to, or strictly greater 
than the critical density k a , c , then we call the traffic state as strictly under-critical (SUC), critical (C), 
or strictly over-critical (SOC), respectively. An under-critical state (UC) can be SUC or C, and an 
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over-critical state (OC) can be SOC or C. Note that different links can have different fundamental 
diagrams. 

From traffic conservation, we have the following system of LWR models in the network [37; 43] 

dK + dQM = aeAyjB (l) 
at dx a 

On a road link, it is well known that (1) admits weak solutions, among which the unique, physical 
solutions should satisfy so-called entropy conditions. For the LWR model, the traditional Lax's 
entropy condition [33] is consistent with vehicles' acceleration and deceleration behaviors [1]. 

On the network shown in Figure 2, (1) is a system of m + n hyperbolic conservation laws, which 
can be understood as a semigroup based on the Godunov method and a Riemann solver [6]. Thus the 
main challenge of network kinematic wave theories is to develop a Riemann solver for the network, 
where hyperbolic conservation equations are coupled with each other at the junction due to merging 
and diverging conflicts among traffic streams. In the Riemann problem, we are interested in finding 
k a {x a ,t) at any (x a ,t) from (1) with the following initial conditions with jump discontinuities at the 
junction: 

k a (x a ,0) = k a , x a <0, aeA (2a) 
k b (x b ,0) = k b , x b >0, beB (2b) 

For upstream link a G A, we introduce a test function <j> a (x a ,t) : (— °°,0] x [0,°°), which is 
smooth with compact support on (— °°,0] x [0,°°); for downstream link b G B, we introduce a test 
function <j) b (x b ,t) : (0,°°] x [0,°°), which is smooth with compact support on (0,°°] x [0,°°). These 

test functions are also smooth across the junction; i.e., (j> a (0,t) = 0&(O,f), and ^(0,f) = ^(0,f) 
for a G A and b G B. Then we attempt to find the weak solutions to the Riemann problem of (1) 
with (2) in the following sense [9]: 

I f I" (k a d ^ + Q a (kAdx a dt+ 



beB 



oo poo 



£ L L (kt^ + QM^dxhdt = 0. (3) 
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2.1 Holden and Risebro's Riemann solver 

In [25], the Riemann problem of (1) with (2) was decoupled into m + n Riemann problems of the 
LWR model based on the following observations. (1) Due to similarity in Riemann solutions; i.e., 
k a {xa,t) = Pa(y), a stationary state, k*, initiates at the junction and spreads on linka (a G A U5). 
That is, for anyx a G (— °°,0] (a G A) orx a G [0,°°) (a G B), lim^oo^ a (x a ,?) = k*, and the stationary 
state pervades the whole link after a long time. Note that it is possible that the stationary state is the 
same as the initial state. It is stationary in the sense that the boundary flux of link a is constant and 
equals 

fa = Qa(K). 
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In addition, due to traffic conservation at the junction, we have 

^ fa = ^ fb-, 

aeA beB 

which is also the Rankine-Hugoniot condition for the junction. (2) A shock or rarefaction wave 
can develop on each link and solves the Riemann problem for the corresponding LWR model with 
initial and stationary states as initial data pairs. That is, for upstream link a (a e A), the shock or 

rarefaction wave solves + d&(*g) = with k a (x a , 0) = < f"' %a < ^ ; for downstream link b 

dt dx " ay a, j | ^ Xfl > Q ' 

(b E B), the shock or rarefaction wave solves ^ + = with k b (x b , 0) = j J*' < ° . 

Thus the Riemann problem is uniquely solved if and only if the stationary states are, and Holden 
and Risebro's Riemann solver is equivalent to finding the following mapping [15] 

(^,...,^ +n ) = R§(^,...,^ m+M ). (4) 

Since stationary states can only propagate backward on upstream links and forward on down- 
stream links, feasible regions of stationary states can be obtained by analyzing the m + n Riemann 
problems of the LWR model. Furthermore, as shown in [25], there can exist multiple feasible 
solutions of stationary states. Thus an additional entropy condition has to be introduced to give a 
unique solution to the Riemann problem of (1) with (2). In [25], an entropy of a junction is defined 
by E = La=" s{^-)^ an d it was proved that there exists a well-defined Riemann solver when g(-) 
is a strictly concave function. However, the entropy condition is not explicitly related to physical 
merging and diverging rules, and vehicles' pre-defined route choices are not considered. 

In [9], vehicles' pre-defined route choices are introduced by a matrix of turning proportions 
Za^b (Vfl EA.be B), where ^ b e [0, 1], and 

£^ = 1, aeA. (5) 

beB 

Therefore, traffic conservation at the junction leads to n equations: Y,aeA fa^a^b = fb Q> e B). It 
was shown that the Riemann solver is well-defined for certain turning proportions. But this Riemann 
solver is not well-defined when m> n. 

In [16; 17; 23; 21; 22], more Riemann solvers along this line have been proposed, but the 
entropy conditions are not directly related to merging and diverging behaviors in these studies. 



2.2 Discrete CTM flux functions 

In [1 1] and [34], the Godunov discrete form of the LWR model was extended to compute traffic 
flows through merging, diverging, and general junctions. In CTM, so-called traffic demand and 
supply, d a (x a ,t) and s a (x a ,t), are defined as functions of traffic density [14; 11; 34] 

d a = D a (k a ) = Qa(rmn{k a ,k a , c }), (6a) 
s a = S a (k a ) = Q a (max{k a ,k ajC }), (6b) 
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where traffic demand increases in total density k a , and traffic supply decreases in total density 
k a . Furthermore, q a = mm{d a ,s a }, and C a = max{J a ,5 a }. In addition, since d a /s a is a strictly 
increasing function of k a , k a can be uniquely determined by d a /s a . We denote this function by 

k a = K a (d a /s a ). (7) 

That is, the pair of demand and supply can uniquely determine the traffic state at a location and 
time. 

Within the framework of CTM, boundary fluxes through a junction at time t are be calculated 
as follows: 

(/l) ' ' - ifm+n) = FF°(<ij,- • • , d' m ,^ m+l , ■ ■ -j^m+n), (8) 

where d[ , • • • , d' m are traffic demands in upstream cells next to the junction, s' m+l , • • • , s' m+n are traffic 
supplies in downstream cells next to the junction, and f[, • • • , f' m+n are the out- and in-fluxes of all 
links. The discrete flux function (8) has been used to model traffic dynamics at various bottlenecks 
within the framework of CTM. For example, when vehicles have pre-defined route choices with 
turning proportions given in (5), a flux function was derived from the First-In-First-Out (FIFO) 
diverging and fair merging rules [32]. Other CTM merging and diverging models can be found in 
[28; 27]. 

Numerically, the flux functions (8) can be incorporated into discrete traffic conservation equa- 
tions and simulate traffic dynamics in a road network with given initial and boundary conditions. 
Thus (8) can be considered as approximate, numerical Riemann solvers [36]. We can see that, using 
the concepts of demand and supply, it is rather straightforward to construct such flux functions for 
solving (1) numerically. While in many other systems of hyperbolic conservation laws, one has 
to solve the Riemann problem first and then obtain boundary fluxes as in the Godunov method. 
Physically, (8) models how the right of way at the junction is allocated among competing traffic 
streams and represents drivers' macroscopic merging and diverging behaviors. Some of these 
macroscopic diverging and merging rules have been verified through observations [38; 35; 40; 8; 3]. 

3 A new Riemann solver 

Since discrete flux functions (8) are both numerically and physically well-defined, it is reasonable 
to use it as an entropy condition to pick out unique, physical solution to the Riemann problem of (1) 
with (2). In [30; 28; 27; 29], a new analytical framework was proposed by replacing the entropy 
condition in [25] by various CTM flux functions in (8). In this section, we first describe the new 
framework and then focus on a new flux function derived in [29]. 

3.1 Riemann solvers with CTM flux functions as entropy conditions 

In order to use CTM flux functions (8) as entropy conditions, we first enlarge the weak solution 
space by introducing a new interior state on each link in the Riemann solutions. That is, in the new 
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weak solutions, we assume thatk a (0,t) = k®, but lim r ^ TC fc a (x a ,?) = k* a for x a G (— °°,0) {a E A) or 
x a E (0,o°) (b E B). Such an interior state is right adjacent to the junction and takes a space of 
measure zero. Thus the new weak solutions to the Riemann problem still satisfy the traditional 
definition in (3). In addition, they do not impact shock or rarefaction waves on all links, which are 
determined by stationary and initial states. 

Such interior states were first observed in numerical solutions of the Burgers equation [47; 7] 
and the LWR model [31] when stationary shock waves occur. Theoretically, if the LWR model 
ff + = on a line x € (— oo ? oo) is solved by a zero-speed shock wave, then an interior state 
can be introduced at the stationary discontinuity without violating the definition of weak solutions or 

( k L , x<0 

the Lax entropy condition. That is, k(x, t) = < kiu, x = with ki<kM < kR and Q(ki) = Q(kR) 

[ k R , x>0 

is a feasible solution with a flimsy interior state &m at x = 0. In reality, such an interior state can 
be observed at the interface of a stationary shock wave, if a detector covers a part of the upstream 
traffic stream and a part of the downstream traffic stream. Depending on the relative location of the 
detector, such an interior state may not be unique. Thus the interior states can physically exist, and 
the enlarged function space of weak solutions is still well defined. In addition, in [28; 27; 29], it was 
shown that the inclusion of interior states is necessary for Riemann solvers to be well-defined with 
some CTM flux functions. Thus the introduction of interior states into weak solutions makes both 
physical and mathematical senses. From the viewpoint of traffic flow modeling, the new framework 
is much more powerful and flexible, since it allows many flux functions derived from driving rules 
as entropy conditions. 

Since discrete flux functions (8) are defined in demands and supplies, it is reasonable to 
use demand and supply, instead of density, as state variables in the new Riemann solver; i.e., 
traffic condition at (x a ,t) is determined by U a = (d a ,s a ). Thus a traffic state is UC if and only if 
d a <Sa = C a , or equivalently U a = {q a ,C a ); a traffic state is UC if and only if s a < d a = C a , or 
equivalently U a = {C a ,q a )- 

In the demand-supply space, we denote the initial and stationary traffic states on link a by (d a ,s a ) 
and (d*,s*), respectively. Then we have the following observations in the Riemann solutions: 

f a = minKV:}, aeAUB (9a) 

In addition, we have the following lemma regarding the feasible regions of stationary states. 

Lemma 3.1 [29] In demand-supply space, the feasible regions of stationary states are given by 

U* a E M(Uar) = (d a ,C a )U{(C a , S * a )\ S * a <da}, aeA (9b) 
U* b E ¥(.,U b ) = (C b , Sb )U{(d* b ,C b )\d* b < Sb }, b E B (9c) 

(9d) 

which lead to 

fa < d a , aeA (9e) 
fb < s b , bEB (9f) 
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Further, we denote the interior state on link a by [/" {a GAUfi). Since interior states do not 
propagate into road links, the Riemann problems with stationary and interior states as initial data are 
solved by waves, whose speeds are positive on the upstream links and negative on the downstream 
links. Therefore, in supply-demand space, the feasible regions of interior states are given by 

U° a G F(I/ fl V) = {(C a X)K < d *a = C a } U \s° a > d* a ,d* a <s* a = Cj, (10a) 

uS e M(;Ut) = {(<r b ,c b )\4 < si = c h }u{(d b , s ° b )\d h > 4,4 < d* b = c b }. (iob) 

where a G A and b G B. When U® = U*, interior states do not exist. 

From the feasible regions of stationary and interior states in (9) and (10), we have the following 
corollary. 

Corollary 3.2 For upstream link a G A, U* is SOC if and only if f a < d a , and U* is UC if and only 
if f a = d a . In addition, ifU* is SOC, then U* = U® = (C a ,f a ); ifU* is UC, then U* = (d a ,C a ), and 
U® — (d®, s®) with s a > d a . For downstream link b G B, U b is SUC if and only if ' f b < s b , and U b 
is OC if and only if f b = s b . In addition, if VI is SUC, then Uj* = U® = (fb,C b ); ifU b is OC, then 
U£ = (C b ,s b ), and Ujj = (d$,$) with d° b > s b . 

Here we define four ratios for upstream link a G A: the initial demand level 8 a = the 

stationary demand level 8* = t ^ L , the interior demand level 8® — ^, and the flux level 8 a = 
Then we have the following results. 

Corollary 3.3 U* for a G A is SOC if and only if 8 a < 8 a ; and U* is UC if and only if 8 a = 8 a . In 
addition, ifU* is SOC, then 

l = S° a =S:>S a >S a ; 

ifU* is UC, then 

K =&a = §a- 

Note that the relationship between 8® and other demand levels when U* is UC is to be determined. 

Similarly we define four ratios for downstream link b G B: the initial supply level o b = g|, the 

s* n s° f 

stationary supply level o b = ^, the interior supply level a b = and the flux level o b = Then 

we have the following results. 

Corollary 3.4 U b for b G B is SUC if and only ifa b < o b ; and U b is OC if and only if o b = o b . In 
addition, ifU£ is SUC, then 

l = o b =a b >o b > 6 b ; 

ifUfr* is OC, then 

a b =o b = d b . 
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Note that the relationship between <j£ and other supply levels when U£ is OC is to be determined. 

In the new Riemann solvers, various CTM flux functions in (8) are used to determine fluxes 
locally from interior states as follows: 

(/l,---,/m+n)=FFVl,---,^,^+l,---,0- (ID 

From the corollaries above, we can see that, given initial upstream demands, d a (a € A), and 
downstream supplies, si, (b G B), if one can find all boundary fluxes, f a and f h , then the stationary 
states can be uniquely determined. But the interior states may not be uniquely determined. Thus the 
Riemann problem is uniquely solved in the sense that all stationary states and waves on all links are 
uniquely solved. Therefore, the Riemann problem is uniquely solved if and only if, from (9), (10), 
and (11), we can find a unique flux function that maps initial conditions into boundary fluxes 

(/lr •• > fm+n) = FF(di, • • • ,d m ,s m+ i, ■ • ■ ,s m+n ). (12) 

Hereafter, we refer to FF (• • •) in (1 1) as local or discrete flux functions and FF( - • •) in (12) as 
global or continuous flux functions. Obviously, FF(- • •) is the Godunov flux function, since it is 
derived from Riemann solutions. 

Then from Corollary 3.2, we obtain the corresponding demand-supply Riemann solver, which 
maps upstream demands and downstream supplies into stationary states: 

(Uf , • • • ,U m+n ) = MS'(d\, ■ ■ ■ , d m ,s m+ i, ■ ■ ■ ,s m+n ). (13) 

Since both demand and supply are many-to-one functions of density, it suggests that different 
initial densities could lead to the same stationary states, if and only if the corresponding upstream 
demands and downstream supplies are the same. From RS'(- • •) we can obtain the traditional Rie- 
mann solver RS(- • •) in (4) by first converting (ki, • • • ,k m+n ) into (d\,-- ■ , d m ,s m+ i, ■ ■ ■ ,s m + n ), then 
applying (13) to find stationary demands and supplies, and finally converting (Z7f , • • • , U m+n ) into 
(k\ , • • • , k* m+n ) using (7). However, from the traditional Riemann solver one may not be able to obtain 
a demand-supply Riemann solver. It can be seen that a Riemann solver in (13) satisfies the consis- 
tency condition [15], if and only if FF(J b • • -,d m ,s m+h - ■ ■ ,s m+n ) = FF(Jj, • • -,d* m ,s* m+l , ■ ■ ■ ,s* m+n ), 
which is equivalent to that (t/*, • • • , U m+n ) = RS'(d*, • • • , d m ,s m+1 , ■ ■ • ,s m+n ). In addition, we call 
the local or discrete flux function FF°(- • •) invariant if FF° (•••)= FF( - • •)• 

3.2 A new local flux function 

In this study, we use the Godunov flux function in [29] as a new local flux function, (11). Note that 
the discrete flux function in [29] is consistent with fair merging and first-in-first-out diverging rules 
but different from the Godunov flux function. Here all vehicles have predefined route choices, and 
the turning proportions, ^ a _j.^, are given in (5). The discrete flux function FF°(- • ■) is defined as 
follows: 
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1. The out-flux of upstream link a E A is given by 

f a = minK,e°C a }, (14a) 
where the interior critical demand level, 0°, is defined as follows: 

6 = mm < max — , mm max >, (14b) 

ya<EAC a beBAiCA LaeAxCa^a^b J 

where Ai is not empty. 

2. The in-flux of downstream link b E B is given by 

fb = E/fl^*- (14c) 

aeA 

In terms of demand and supply levels, FF°(- ■ •) in (14) can be re-written as 

~8 a = min{5 fl °,0 }, (15a) 



„0 sO Cb&b EaeA\Ai Ca^b^a 
6 = mm < max o a , mm max — 

aeA beBAiCA LaeA^Ca^b 



(15b) 



Cb&b = £^4, (15c) 

aeA 

where C a -+b = C a £, a -^b for a G A and b E B. 

3.3 Average demand levels 

In this subsection, we discuss properties of 0° in (15b). For a E A and b E B, we define the 
demand and supply levels by \i a E [0, 1] and E [0, 1], respectively. Further, we denote %b — 
CbVb ~ L«eA Ca^bHa- 1 We denote the vector of \x a for a E A by fl, and the vector of v b forbEB 
by V. We define the average demand level of set A\ for link b by 

N Q>Vft - LaeAVAi C a ^bVa Kb+HaeA^a^bHa .... 

Yb{Ai) = = , (16) 

LaeAi ^a^b LaeAi ^a^rb 

where A\ C A is non-empty. 

Lemma 3.5 For a EA\, andA 2 =A\ \{a} ^ 0, if [i a >,=, < Yb(M), then Y b (A 2 ) <,=,> Yb{M), 
and jl a >, =, < Yb{A 2 ), respectively. For a ^ A\, and A 2 = A\ U {a}, <, =, > %(Ai), 
7^(A2) <,=,> 7^(Ai), and \x a <,=,> 7fo(A2), respectively. That is, if we remove a link with a 
larger demand level, then the average demand level decreases; if we add a link with a larger 
demand level, then the average demand level increases. 



'in this study we do not consider situations when t, a ^b = 0- 
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Proof. For a G A\, and A 2 = Ai \ {a} ^ 0, we have 

, A \ LaeAi Ca^bHa + C a ^bHa + 

7b{A\) = ^ ~ — ; , 

LaeA 2 ^a^b ' ^a^b 

which leads to 

7*( A 2j = = ~ = %( A U H ^ t; » 

LaeA 2 ^a-+b La<=A 2 ^a^b 

Thus, if n a >,=,<y b (A 1 ), then y fc (A 2 ) <,=,> Yb{M), and >,=,< %(A 2 ), respectively. 
For a ^ Ai, and A 2 = A\ U {a}, we have 

/ . \ HaeAy Ca^bHa — C a ^ b Ha + 7t b 

7b(Ai) = 2 - — , 

LaeAi ^a^b ^a^-b 

which leads to 

v La€A 2 C a ^bHa + Kb , {jb{A\) - Ha)C a ^b 

Yb{A2) = = = Yb{Ai) . 

LaeA 2 ^a-^tb La<EA 2 ^a^b 

Thus, if ji a <,=,> Jb{A\), 7&(A 2 ) <,=,> Yb{A\), respectively. Furthermore, 

{YbiM) - Ha) Ca ^ b = H C a^bHa + C a ^bHa + K b - fl a £ C a ^ b -C a ^ b \i a 

aeA 2 a€A\ aeA\ 

= C ^bHa + X b -fla Y C ^b = {Yb{A\) ~ Ha) £ C a ^ b . 

a£A\ flGAi a^zAi 

Thus, if Ha <,=,> Yb(Ai), Ha <,=,> Yb{Ai), respectively. ■ 
We denote the maximum average demand level for link bbyT b = max Al c AAl7 ^) Yb{A\), which 
has the following properties. 

Lemma 3.6 Yb(A\) and T b have the following properties. 

1. When % b > 0, then y b {{a\) > Ha far a E A, and T b = max aeA Yb({ a }) > m & x aeAHa- 

2. When % b = 0, then Yb({ a }) = Ha far a G A, and T b = max aeA Yb({ a }) = m ax aeA Ha- 

3. When 7t b < 0, then Yb({ a }) < Haf° r a & A, and there exists a unique A\ 7^ 0, such that 
Tb = Yb(M), Ha>T b > Ha far a G Ai and aeA\A h 

Proof. 

1. When % b > 0, from (16), we can have Yb({a}) > Ha for aeA. Then from Lemma 3.5 we have 
min fleAl Ha < Yb(Ai) < max aeAl %({a}). Thus T b = max fleA %({a}) > max aeA jU a . Note that 
max aeA Yb({a}) and max aeA Ha may attain their maxima for different a. 
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2. When K b = 0, from (16), we can have Yb{{a}) = Ha for a G A. Then from Lemma 3.5 we 
have min aeAl Ha<Yb(^l)< max aeA] \i a . Thus T b = max aeA %{{a}) = max fleAi u fl . 

3. When Kb < 0, from (16), we can have Yb({a}) < Ha f° r aeA. Since A has a finite number 
of subsets A\,Tb exists and we can find A\, such that Y b = 7,(Ai), and the value of T b 
is unique. First, for a A\, Ha < Yb{A\), since, otherwise, from Lemma 3.5 we have 
Yb{A\) < Yb{A\ U {cc}), which contradicts that = 7,(Ai). Second, there exists at least one 
a G A\, such that \x a > Yb{A\), since, otherwise, from Lemma 3.5 we have Yb({a}) = Ha 
for all a G A\, which contradicts that Yb({a}) < Ha- Third, for all a G A\, \l a > Yb{A\), 
since, otherwise, from Lemma 3.5 we have Yb{A\) < Yb{A\ \ {a}), which contradicts that 
Tfc = Yb(Ai). Assume that a G A\ c A\ and fi a > Yb{M), then ji a = Yb{M) for a G A\ \A\, 
and Yb{A\) = Yb{A\)- Since A\ ^ is unique, there exists a unique A\ = A\ ^ 0, such that 

= 7>( A i)> M« > r & > M« for a G Ai and a gA\Al 

■ 

Assuming that \i a (a G A) are in a decreasing order; i.e., Hi > • • • > \l m , we define the following 
average demand level of the first / upstream links: 

Yb(l) = Yb({h -,/})■ 
Then we have the following corollary. 
Corollary 3.7 %(/) and Tb have the following properties: 

1. When K b >0, Yb(l) > Hi far I e A, 7^(1) > ••• > Yb(m), and Y b = max aeA Yb({a}) > 7>(!) > 
Mi. 

2. Wzen % = 0, > M// or f £ A, %(1) > • • • > 7,(ra), and T fo = Hi- 

3. When K b < 0, there exists a unique I* G {1, • • • , m}, suc/i ?/za? 7,(1) < • • • < Ybi)*) > 7>(^* + 
1) > ••■ > 7>(w), 7,(Z) < M//or/= 1, •••,/*, andYb{l) > Hifa r I = ^* + V ' • ? m - That is, for 
every b, there exists a unique 11 solving min/, smc/i ?/za? 7,(7) > 7,(/ + 1). In addition, l£ = 1 
and 7,(1) > Ml if and only if Kb > 0. 

Proof When K b > 0, the results are obvious. Here we will focus on Kb < 0. We denote I* as the 
number of links in A\, where Y b = 7,(Ai) and Ha > Tb > Ma f° r a G A\ and a G A \A\. Then /* 
is unique, and Hi* >^b = Yb(l*) > Hb(l* + !)• From Lemma 3.5, we have 7,(1) < • • • < 7,(/*) > 
%(/* + l) > • > %(m), 7,(/) < M/ for / = 1, •••,/*, and %(/) > jU Z for / = /* + 1, • • • ,ra. ■ 
We define the critical demand level by 

/ N . CbVb — LaeA\Ai C a ^bHa 
g(jl,V) = minmax ^ 

b A l LaeAi ^a-^b 

K b + LaeA, C a ^bHa , 1 

= minmax (17) 

b A l LaeAi ^a^b 

Since v) = min^r^, we have the following theorem regarding g(fl 7 V). 
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Theorem 3.8 For C a ^ b >0,C b > 0, \x a G [0, 1], and v b G [0, 1] (a G A and b G B), g(fi, v) defined 
in (17) has the following properties: 

1. g(fl, v) > max fl \i a if and only if min b<E B n b > 0. 

2. If and only if min be B7t b < 0, there exists a unique non-empty A* C A such that 

g(fl,V) = minmax%(Ai) = min%(A*) = maxmin7,(Ai), (18) 

b A\ b A\ b 

and 

min^ > g(jl,v) > max \x a . 

aeA* aeA\A* 

Proof. When min^ 7T& > 0, from Lemma 3.6 we have v) = min^ T b > max a \i a . 

When min^ 7Z^ < 0, we denote A b = {a E A\/J, a > T b }. From Lemma 3.6, A ft = if and only if 
7t b > 0. We denote A* = U^ g bA^ ^ 0, since min^ ^ < 0. In addition, for any a G A*, /i a > min^ T^, = 
g(|i,v); and for any a G A \A*, jU a < min^r^, = g(|i,v). In addition, g(fl,v) = min b y b (A*), and 
min fleA , jU a > min fc %(A*) > max aeA \ A ^ \i a . 

Since min b n b < 0, we have mm b y b ({a}) < fl a for a G A. Since A has a finite number of 
subsets Ai, thus we can find A\ such that max Al rmn b y b (Ai) = mm b y b (A\), and the maximum 
value is unique. First, for any a £ A\ fl a < rcm\ b y b (A\), since, otherwise, from Lemma 3.5 
mm b y b (A\ U {a}) > mm b y b (A\). Thus a G A\ if fi a > mim, y fe (Aj). Second, there exists at least 
one a G A*, such that /i a > min^ 7&(Aj), since, otherwise, jU a < min^ y b {A\) for all aEA\, and from 
Lemma 3.5 fl a < mm b y b ({a}), which is not possible. Third, if a G A\, then /i a > min^^A*), 
since, otherwise, from Lemma 3.5 mm b y b (A\ \ {a}) > mm b y b (A\). Therefore, without loss of 
generality, we can remove all a from A\ if \i a = min b y b (A\), and the value of min^ y b (A\) does not 
change. Then A* = {a G A\fi a > min b y b (A\)} is unique; that is, there exists a unique A\ such that 
min aeA *jU a > mm b y b (A\) > max aeA \ A .ju a . 

Therefore, A* = A\, and (18) is proved. ■ 

From Theorem 3.8 we have the following corollary. 

Corollary 3.9 If jl a is decreasingly ordered, we have the following results on g(fl,v): 

1. If and only if min b n b > 0, gQll, V) > jUi. 

2. If and only if mm b 7t b < 0, there exists a unique I* G {1 , • • • , m}, such that 

g(H,v) = min max 7,(7) = min7,(/*) = max min 7,(7), (19) 

b l b lb 

and 

Hi* >g(p,v) >jU/*+i- (20) 
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Figure 3: The pattern of Yb{a) 



We show an example in Figure 3, where jUi > • • • > /iio- For b = 1, 7i(l) < ■ • ■ < 7i(6) > 
7l(7) > ••• > 7i(10), 72(a) < \L a for a = l,---,6, and 72(a) > jU a for a = 7, •••,10. For & = 2, 
72(1) < 7 2 (2) = 72(3) = 7 2 (4) > ••• > 72(10), 72(a) < \x a for a = 1,2, and 72(a) > \l a for a = 
3, • • • , 10. For b = 3 75(1) > • • ■ > 73(10), and 73(a) > \l a for a = 1, • • • , 10. Then we can verify the 
lemmas above and find that /* = 6, such that /i6 > g(fi, v) = 71 (6) > ^7. 

Since 0° = min{max aeA 8°,g(8°, <T )}, where 5° = (8°) aeA and <T° = (a°) foefi , we have the 
following observations from Theorem 3.8: 

1 . When C h ol > £ aeA C a ^ h 8° for all b E B, then g(S°, a ) > max aeA 5°, and 0° = max fleA 5°. 

2. When 0,0^ < LaelCa-^a f° r some b E B, then there exists a unique, non-empty A* such 
that 0° = g(8 ,O°) = min b y b (A*), and min fleA „ <5 fl ° > 0° > max ael \ At 5°. 

Thus 0° is well-defined, bounded between and max aeA <5°, and continuous in (5°, ■ • ■ , 5^, Oj^ +1 , ■ ■ • , o^ +n ) • 



4 Solutions to the Riemann problem 

From the preceding section, we can see that the demand and supply levels satisfy the following 
conditions: 

{1 = 5° = 5; > 5 fl > 4} or = 8 a = 5 a and s° > 4} (21a) 
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{ 1 = ctJ = ol > e b > d b } or {<7* = o b = a b and s° b >d b } (21b) 

5 fl = min{5°,0 } (21c) 

0° = max aeA 5 fl ° or (5 fl °) fleA ,^ > 0° = min, efi $(A*) > (5°) fl€AV4 , (21d) 

Q,Ofc = E aeA (21e) 

where 7? (A*) = c ^ C «^/A ^ j n this section, we attempt to solve S a (a e A) in 8 a and a b ; 

i.e., a mapping from upstream demand levels and downstream supply levels to upstream flux levels: 

(<5i,---,<5 ra ) = FF / (5i,---,5 ra ,a m+ i,---,a m+ „). (22) 

Then from the definition of 8 a and (14c) we can have the flux function (/i , • • • , f m +n) = FF(di , • ■ ■ , J ra , s m+ i , • • • , s m + M ) 
and the corresponding Riemann solver. 

4.1 Further properties of demand and supply levels 

From Corollary 3.3 and (21) we have the following lemma. 

Lemma 4.1 U* is SOC if and only if 6° = 8 a < 8 a , and U* is UCifand only if 6 > 8 a - In addition, 
if 0° < 8 a , then U* = U% = (C a , f a ) is SOC, and 

\ = 8 a ) = 8:>8 a >8 a = e Q - 

ifO° > 8 a , then U* = U% = {d a ,C a ) is SUC, and 

8° = 8* a =8 a = 8 a <6°; 
if 8° = 8 a , then U* = (d a ,C a ) is UC, = (d^s° a ) with d Q a > d a and s° a > d a , and 

8°>8: = 8a = 8 a = d°. 

In the third case, U® = U* if and only if 8® = 0°, and the interior state U® is different from the 
stationary state U* when 8® > 6°. 

An example of the relationships between 8 a , 8*, 8®, and 0° is shown in Figure 4, in which U* = U® 
is SOC for a = 1, • • • ,4, U* = U° is UC for a = 6, • • • , 10, and U* is UC but U* ^ £/ 5 °. 

We define A° = {a e A|5 fl ° > 0°}. From Lemma 4.1, U* is either SOC or UC with U* ^ t/ fl ° for 
a G A®; and d^ = d a for a e A \A°. From Theorem 3.8, we can see that min^ eB 71® = mm be BS b — 
l% =1 C a ^ b 8% > iff AO is empty. If A° is non-empty, then (5°) a6A o > 0° = min, eB ^(A°) > 
(<5„) aeA \ A o. Then we have the following lemma. 

Lemma 4.2 For downstream link b, ifA®^% and ^(A°J = 0° or if A® = and 7t® = 0, then 
I/]* = U® = (C b ,s b ) is OC, and 

°b = °b =°b = & b \ 
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Figure 4: The relationships between 8 a , 8*, 8®, 8 a , and 0° 



When A° ^ and $(A°) > 0° or when A° = and 7$ > 0, there can be two types of solutions: (i) 
U* = Ujj = (f b ,C b ) is SUC, and 

l = o$ = o* b >o b > 6 b ; 

(ii) if f b = s b < C b , another possible solution is that U£ = (Q,, s b ) is SOC,U^ = (dls° b )^U;with 
d? > s b and s9 > s b , and 

o£ > o£ = o b = a b . 

Therefore, ifU* is SUC, A° = and > or A° ^ and $(A°) > 0°; ifU* is OC, the interior 
state may not be the same as the stationary state. In addition, for all b, 6 b < o b < (7 b < o b , where 
all the equality signs hold when A® = and s® = EaeA C a ^ b 8^, or when A® ^ and 0° = y^(A®). 

Proof. 

1. When A° ^ and ^(A°) = 0°, then from (14c) we have 

f b = 0° £ C a ^ b + £ C a ^ b 8°(a) 
aeA? aeA\A° 

- ^ L C a-^b+ L L a^bO {<X)—S b , 

LaeA° ^a^b aeA o a<EA\A° 

If f b < s b ; i.e., U b =U b = {f b , C b ) is SUC, and from (14c) we have f b = s Q b = C b < s b , which 
is impossible. Thus f b = s b = s b , and U b = U b — (C b ,s b ) is OC. 
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2. When A° = and 7t b = 0, then from (14c) we have f b = s b . Similarly, we have fb = Sb = s b , 
and [/* = £/° = (C b ,s b ) is OC. 

3. When A° ^ and $(A°) > 0°, then from (14c) we have f b < s° b . If f b < s b ; i.e., = Ujj = 
(f b ,C b ) is SUC, then s° b = s* b = C b > s b > f b , and 1 = o° = >a b > c b . If f b = s b ; i.e., 
U b = (C b , s b ) is OC, and Ufj = (d^ s^) ^ U b with d b > s b and s b > s b , which is possible only 
if s b < C b . In this case, o b > o b = o b = d b . 

4. When A° = and n b > 0, similarly, there are two types of solutions: U b = U b = (f b ,C b ) 
is SUC, then s° b = s* b = C b > s b > f b , and 1 = oj = o£ > a b > a b ; = (C b ,s b ) is OC, 
U b = (d b ,s b ) 7^ U b with d b > s b and s b > s b , which is possible only if s b < C b , and o b > 
° b = a b = o b . 



4.2 Solutions of the flux function 

In addition to the interior critical demand function 0° in (14b), we also define three other critical 
demand functions as follows: 

= minjmax o a ,min max — j—z } = minjmax o a ,g(o,<7)}, 

aeA b A\ LaeAiCaha^b aeA 

= minjmax 8 a , mm max — } = minjmax 5 a ,g(8 ,G )}, 

aeA b Ai LaeAi^a^a^b aeA 

Z . r ~ s b ~ HaeAXA, fa^a^b , . r =; ,x 

= minjmax o a ,min max —z }= minjmax o a ,g{o,0)$. 

aeA b A\ LaeAiCaha^b aeA 

We define the residue supply of link b by 7t b = s b — E™=i d a ^ a -^b- 

In the following, we demonstrate that there exists a global flux function, (12), satisfying (9), 
(10), and (14). We also show that the local flux function in (14) is invariant. 

Theorem 4.3 For the Riemann problem, we have 

0° = = e* = O, (23) 
8 a = mm{8 a ,e} = rmn{8? l ,e°} = mm{8:,e*} = rmn{8 a ,e}. (24) 

Therefore, the discrete flux function (14) is invariant. 

Proof First, if min b K b > 0, all upstream links are stationary at UC; i.e., q a = d a for aeA. Otherwise, 
from Lemma 4.1 we have that 0° < max a 8 a < max fl <5^, and A* = ja|<5^ > 0°} is not empty. From 
Theorem 3.8, there exists b, such that )^(A*) = 0°. Further from Lemma 4.2 we have that f b = s b . 
However, from (14c) we have the following contradiction: 

fb = d ° C a+ d a^a^b<Y, da ^b<Sb, 

aeA* aeA\A* aeA 
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since 0° < max a 8 a and % b > 0. Thus, no upstream links can have SOC stationary states, and 
from Lemma 4.1 we have that either 8® = 8* = 8 a = 8 a < 0° or 8® > 8* = 8 a = 8 a = 6° for 
any a E A. Thus when A* = {a\S® > 0°} is not empty, then 0° = max aeA <5 fl = max aeA <5* — 
max aeA 8 a < max aeA <5° 2 ; otherwise, 0° = max aeA 8 a = max aeA 8* = max aeA 4 = max aeA <5°. 3 
In addition, since s* h > s b > f b , s* b - YZ=\ d ^a-^b > 0, and s b - YZ=\fc&a-*b- From Theorem 
3.8, = max aeA <5 a , 0* = max aeA <5^, and = max fleA 5 a . Therefore, = 0* = 0° = 0, and 
8 a = min{S a ,0} = min{5;,0*} = min{5 fl °,0 } = min{5 fl ,0}. 

Second, if min b it b < 0, then at least one upstream has a SOC stationary state. Otherwise, f a = d a 
for all a, from (14c) we have that f b = Y^aeA^a^a^b > s b when n b < 0. This contradicts f b < s b 
for all b. From Lemma 4.1, A® = {a\0° < 8®} is not empty. Thus 8 a = 0° for a E A*, and <5„ = 
8 a = 8a = 8 a for a E A \A®. From Theorem 3.8, we have min aeA o 8® > 9° = mm beB y^(A°) > 

max aeA\A0 5 a> where 

M) 



s b ~ £ceeA\A° C a ^ b 8 a _ C b O b ~ EaeA\A° C a ^ b 8 a 
EaeA° Ca->b EaeA° C a ^b 



In addition, 



, N Cb ° b ~ ^aeA\A« C a ^ b 8 a 

7b\A^) — ~ — s 

LaeA® ^a-^b 

^tf A 0\ _ s b~ £aeA\A° C a ^b8 a ^ C b O b — LaeA\A« C a ^ b 8 a 

2-aeA2 ^a^b LaeA° ^a^b 

s b ~ LaeA\A$ C a ^ b 8 a C b O b — L«ga\ao C a ^ b 8 a 



From Lemma 4.2, there exists b such that y$(A®) = 0°, for which o b = o b = o b = o b , and 0° = 
YbM = ft(A°) = ^(A»); for other b E B, we have that $(A«) > 0°, and f b = e°£ aeA oC a ^ b + 
LaeA\A»C a ^b8a < s b < s* b , which leads to 0° < 7 fc (A«) = 7 fc (A«) < ^(A«). Thus we have 

min y b (A°) = min7^(A°) = min y b (A°) = min^(A^) = 0°. 
beB ' v *' beB lbK * ; beB 1 K *' beB lb * ; 

If we denote A* = {a\6° < 8 a }, which is the set of upstream links with SOC stationary 
states, then A* is not empty. From Lemma 4.1, we have that A* C A^ and 0° = 8 a = 8 a = 
8* < 8® for a E A® \A*, which is the set of all upstream links with interior states. From 
Lemma 4.1, y b {A*) = f b (A,) = 0° when %(A°) = f b (A°) = 0°; and y b (A*) = f b {A*) when 
Yb(A*) = YliA®) > 0°. Therefore, mm beB y b (A*) = mm b(EB yl(A*) = 0°. From Lemma 4.1 we 
have that 0° < 8 a < 8* for aEA*, and 0° > 8 a — 8* for a E A \ A*. Thus we have mm aeAf 8* > 



2 In this case, there exists at least an interior state which is different from the corresponding stationary state. 
3 It is impossible that 0° > ma\ a£A 5„ from the definition of 0°. 
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min aeAt 8 a > mm beB %(A*) = min beB f b (A*) = 0° > max aeA \ A ^ 8 a = max aeA \ A ^ 8*, which leads 
to min beB y b (A*) = = 6* = 0° from Theorem 3.8. Since s b >f b = £™ =1 C a ^ b 8 a , we have from 
Theorem 3.8 that = max aeA 8 a = 0° since A^ is non-empty. In addition, when a E A \A®, 
0° > 8 a = 8 a = 8 a = a ; when a e A°, 0° < 8 a < 8* < 5°. In both cases, (24) is satisfied. ■ 
Remark. Note that Theorem 4.3 is also true when we switch the min and max operators in 
definitions of 0, 0°, 0*, and 0. From the proof of Theorem 3.8, we can see that these values 
remain the same when mm b % b < 0. But if min/, % b > 0, the min- max and max-min operators could 
yield different values, which are not smaller than max a 8(a). A counter example is as follows: 



d= [0.5,0.5], C= [1,1], % 



0.2 
0.2 



then 7 : 



0.75 1.5 0.7 
1.5 0.75 0.7 



In this 



0.8 0.2 
0.2 0.8 

case, min^maxAj y b {A\) = 1.5, and max^ rmn b y b {A\) = 0.75. But both values are greater than 0.5, 
the maximum demand level. Thus Theorem 4.3 still holds. 

Remark. From Theorem 4.3, we can find a unique flux function (f\ , • • • , f m + n ) = FF( d\ , • • • , d m ,s m +\ , • • • , s m+n ) 
defined in (12) and solve the Riemann problem in the following steps: 

1. Calculate from initial conditions in d a , s b , and ^ a ^b- 

2. Calculate f a = min{d a , 0C a }, and f b = LaeAfa^b- 

3. Determine stationary states and interior states 4 from Corollary 3.2. 

If we denote = g(5, a), then 0=0 when mm blEB n b < 0, and > when min beB n b > 0. 
We can see that f a = min{<i a ,0C a }. We define s+ = &C a . Then q a = min {d a ,s~^}, and s^ can 
be considered as effective downstream supply of upstream link a. We define & b by Q_ b = 
rmn J 6es\{i} max AiCA7j3(Ai) 5 . Obviously ® b > 0. We have the following lemma. 

Lemma 4.4 The in-flux of downstream link b can be written as 

fb = mm{dj;,s b }, 

where 

d b = E min H> -* c «}<^- 

aeA 

That is, d b can be considered as effective upstream demand of downstream link b. 

Proof. When min beB n b > 0, from Theorem 3.8 we have that > max aeA 8 a , and & b > max aeA 8 a . 
Thus d h = YsaeAda^a^b < $b since Kb > 0. From Theorem 4.3 we have f b = d b = min{d b ,s b }. 

When mm beB n b < 0, if & b = 0, then f b = d b = £aeAfa£a->b < s b\ if ®-b > ®> then from 
Theorem 3.8 there exists A* ^ such that yb(A*) = and mm a<EA ^ 8 a > > max aeA ^ s 8 a . There- 
fore, 

Sb = ®J^Ca£a-+b+ E d a^b=Y min i d ^ @C aKa^b = fb- 

a<EA* aeA\A* aeA 



4 The interior states may not be unique. 

5 Note that if b is the only downstream link, then d-t is set to 1 . 
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In addition, we have 

fb = min H' ® C a}Za-+b < Y min H> ®-bC a }£a^b = d h . 

Therefore, = min{<i^ , Sf,}. ■ 
From Theorem 4.3 and Lemma 4.1, we have the following corollary. 

Corollary 4.5 U* is SOC if and only if d a < s+, and U* is UC if and only if d a > s£. In addition, 
ifd a > Sq, then U* = U® = (C a ,q a ) is SOC, and 

l = 8° = 8Z>8 a >8 a = e; 

if d a < S+, then U* = U® = (d a ,C a ) is UC, and 

8 Q a =8: = 8 a = 8 a <d; 
ifd a = s t' tnen Ua = (d a ,C a ) is UC, U® = (d®,s®) with d® > d a and s® > d a , and 

8° a >8: = 8 a = 8 a = e. 

In the third case, it is possible that U® ^ U*, but U® = U* is also a valid solution when 8® = 8 a . 

Similarly, from Theorem 4.3, Lemma 4.2, and Lemma 4.4 we have the following corollary 
regarding downstream links. 

Corollary 4.6 U b is SUC if and only if Sb > d b , and U£ is OC if and only if Sb < d b . In addition, 
if Sb > d b , then U b = U b = (q b ,Cb) is SUC, and 

l = o$ = o b *>o b > o h \ 

if Sb < d b , then U b — U® = {sb,Cb) is OC, and 

a b = °b = a b = o b \ 
if Sb = d b , then U b = (sb,Cb) is OC, U b = (s^s®) with s b > s b and s b > Sb, and 

Ob > o£ = o b = o b . 

In the third case, is is possible that U b ^ U b , but U? = U b is also a valid solution when o b = Ob- 
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5 Discussions 



5.1 Special cases 

For a linear junction with m = n = 1, ^x^-2 = 1, and /i = fa. In this case, 0=^,0 = min{<5i, 
and /i = min{0, 5i} = min{0, Thus the effective supply of link 1 is s\ = &C\ = s 2 , and the 
effective demand of link 2 is d 2 = mm{d\ , & 2C1 } =d\. From Theorem 4.3 we have fi = fa = 
mm{dx,s 2 }. Thus from Corollaries 4.5 and 4.6 we have the following three cases. 

1 . When s 2 >d l ,d 1 <s+, and s 2 > d 2 . Thus U{ = t/f = (<*i , Ci ) is UC, and C/| = U$ = (d h C 2 ) 
is SUC. In this case, there is no interior state on either link. 

2. When s 2 < d u d\ > s+, and s 2 < d~ . Thus £/* = £7° = (C h s 2 ) is SOC, and C/| = C/ 2 ° = 
(^2,^2) is OC. In this case, there is no interior state on either link. 

3. When s 2 = d\, d\ = s|, and s 2 = d~. Thus, U{ = (d\,Ci) is UC, and C/| = (C 2 ,s 2 ) is OC. In 
this case, there can be interior states on both links: U® = (d® 7 s®), where d® >d\, and > d\\ 
U 2 = (d 2 ,s 2 ), where d 2 > s 2 , and s 2 > s 2 . However, since fa = minjjp^j}, we have either 
d® = d\ or s 2 = s 2 . Thus Uf = C/f or U 2 = U 2 ; i.e., there can be only one interior state. 

For a merging junction with m > 1 and n = 1, £, a -^m+i = 1, and 

s m+\ — LaeA\Ai 

= max — — . 

A iCA LaeAxCa 

Thus the effective downstream supply for link a e A is 5+ = 0C a , and the effective upstream 
demand for link m+ 1 is J m+1 = LaeA<4> since ^ = 1. Then we can solve for stationary and 
interior states by following Corollaries 4.5 and 4.6. In particular, when m = 2, we have 

f s 3 -d 2 S3-d\ 53 

1 Ci C 2 'Ci+Cz 1 ' 
q x = min{Ji,0Ci}, 

q 2 = min{J 2 ,0C2}, 
which is consistent with the fair merge model in [28]: 

qx = mm{di,max{s 3 -d 2 1 53}}, 

Ci +C 2 

c 2 

q 2 = min{J 2 ,max{5 3 - d\, - - s 3 }}. 

Ci +L 2 

For a diverging junction with m = 1 and n > 1, 

= mm — g . 
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Thus the effective downstream supply for link 1 is = 0Ci = min]^ ^~ h ^ an< ^ tne effective 
upstream demand for link b E B is d b = min{di , min^l^ ^ b c ^ ^ j^i^. Then we can solve for 
stationary and interior states by following Corollaries 4.5 and 4.6. In particular, 



1+M Sb 

q\ = min{<ii , ©Ci } = min{<ii , min }, 

which is consistent with the FIFO diverge model [27]. 
For a junction with m = 2 and n = 2, we have 

= mm{max{— , — — — }, 

<-lS1^4 <-2S2^4 <-lS1^4 +<-2S2^4 

S4 — ^2^4 S4 — d\t,\^4 S4 

m^X-j^ — ^ — ? 5 — ~p, — ? 1 — ? ~p, — ? / /, 

MS1^4 C2S2^4 MS1^4 +C2S2^4 

fx = min{Ji,0Ci}, 
f 2 = min{J 2 ,0C 2 }. 

From Theorem 3.8, we can see that there are the following scenarios: 

1. Both links 1 and 2 are stationary at UC if and only if > max{^-, ^}. 

2. Link 1 is stationary at SOC and link 2 at UC if and only if ^ > > ||. 

3. Both links 1 and 2 are stationary at SOC if and only if < max{^, ^}. 
We can find all stationary and interior states by following Corollaries 4.5 and 4.6. 



5.2 A simplified framework without interior states 

Within the framework defined in Section 3.1 and a flux function in interior states in (14), the 
Riemann solver is well-defined since fluxes can be calculated from Theorem 4.3, and stationary 
and interior states can be determined from Corollaries 4.5 and 4.6. 

From Corollaries 4.5 and 4.6, interior states can be the same as stationary states in all scenarios, 
and we can introduce a simplified framework as follows: 

1. Stationary states arise near the junction on all links and satisfy (9). 

2. An entropy condition is defined by the flux function in stationary states: 

e* = min{maxS;,s(r,<7*)}, 

aeA 

f a = minK,0*C a },VaeA 

Sb = Y,Sa^VbeB 

aeA 
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From Theorem 4.3, the Riemann problem is uniquely solved and the out-flux of upstream link a is 
still given by f a = min{<i a , 6C a }. Further, Corollaries 4.5 and 4.6 can be used to determine unique 
stationary states. 

Note that, however, not all CTM flux functions can be used in the simplified framework. For 
example, with the following flux function proposed in [32] 

d° 

f a = minK°,min- fa—fih ( 2 ^) 

fb = Y^f^b, (25b) 

it was shown that f a = min{J a , 0C a }, but it is possible that interior states are different from 
stationary states. In this case, FF°(j| ) , • • • , J° , s Q m +1 , • • • , s% +n ) = FF(d\ ,---,d m , s m+ i, • • • , s m+n ) = 

FF*(4,--.,j*, s ; +1 ,.--,4 +M ),but 

FF (d l , • • • , d^, s Q m+l , s Q m +n ) ^ FF(d l , • • • , J° , s Q m +1 , • • • , s m+n ) . 

That is, FF (• • •) ^ FF(- • •), and it is not invariant. Therefore, with an invariant flux function, 
the network kinematic wave model (1) can be defined in the function space of traditional weak 
solutions. But for a non-invariant flux function, e.g., (25), the function space has to be extended to 
include interior states. 
Clearly we have that 

RS(MS(fc 1 , • • • , k m+n )) = RSih , • • • , k m+n ). 

That is, the Riemann solver is consistent in the sense of [15]. In addition, all Godunov flux functions 
can be used as entropy conditions in the simplified framework, in which interior states are the same 
as stationary states. 



5.3 A stationary junction network 

A junction network is stationary if and only if all initial states are the same as stationary states; i.e., 
U a = U* for aeA, and U b = U£ for beB. From Corollaries 4.5 and 4.6, it is possible that d a = d* 
(a E A) and s fc = s* h (b G B). Therefore, such stationary states always exist. In addition, if U a and Ub 
are stationary, then 

f a = mm{d a ,®C a } =mm{d a ,s a }, 
f h = mm{d b ,s b } = mm{db,s b }. 

From Corollary 4.5 we can see that, if U a = (C a , s a ) is SOC, then &C a = s a = q a ; if U a = (d a , C a ) 
is UC, then © > 8 a . From Corollary 4.6 we can see that, if Ub = (db,Cb) is SUC, d b = qb < Cb\ if 
U b = (Cb,s b ) is OC, then d b > s b . 
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When a junction network is stationary, we can replace any SOC state (C a , 9C a ) by the corre- 
sponding SUC state (6C ai C a ), and 6 remains the same; but we may not replace a SUC state by 
the corresponding SOC state. Thus for the same fluxes, there can be multiple combinations of 
stationary states. Furthermore, if we replace d a by q a = 0C a for a E A*, then all upstream links will 
be stationary at UC, and the critical demand level and all upstream flow-rates remain unchanged. In 
addition, if all upstream links are stationary at UC, we can replace an OC stationary state on link b, 
(C h ,q b ), by an UC state, (q b ,C b ). 

6 Conclusion 

The kinematic wave model of network traffic flow, (1), is thus well-defined by the following rules. 
Rl. The constitutional law: q a {x a ,t) = k a (x a ,t)v a (x a ,t) at any pointy on link a and time t. 

R2. The fundamental diagram of speed-density relation, v a (x a ,t) = V a (k(x a ,t)), flow-density 
relation, q a (x a ,t) = Q a (k(x a ,t)) = k a (x a ,t)V a (k a (x a ,t)), demand-density relation, d a (x a ,t) = 
Q a (min{k ajC ,k a (x a ,t)}), and supply-density relation, s a (x a ,t) = Q a (max{k ajC ,k a (x a ,t)}). 

R3. Traffic conservation: + §J = 0. 

R4. Weak solutions with interior states: the kinematic wave model can have discontinuous weak 
solutions and interior states at stationary discontinuities. 

R5. Entropy conditions with local flux functions at any junction at x: we denote the set of upstream 
demands by d(x~ ,/) and the set of downstream supplies by s(x + ,t), then the set of boundary 
fluxes f (*, = FF° (d(jc - , t) , s(x+ , t) ) . 

Thus, if initial conditions in densities and boundary conditions in demands at origins and supplies at 
destinations are given, (1) can be uniquely solved with the aforementioned five rules. 

This modeling framework makes the entropy condition explicit by using boundary flux functions. 
In a sense, it is the reverse process of Godunov method, in which flux functions were derived by 
solving Riemann problems with entropy conditions defined in characteristics or other approaches 
[33]. It is possible to extend this framework for more complicated situations for multi-class traffic 
on multi-lane roads or for other types of intersections. The challenges will be related to identifying 
fundamental diagrams and developing boundary flux functions. Therefore, the Riemann solver can 
be used to determine whether a flux function is well-defined both mathematically and physically. 

In addition, it is possible to extend this framework to study other systems of hyperbolic 
conservation laws, in which demand and supply functions are well-defined. For example, numerical 
E-0 flux function and other approximate Riemann solvers can be used as entropy conditions when 
solving the Burgers equation. 

In the future, we will be interested in analyzing traffic dynamics in a road network with the help 
of the new kinematic wave model. Such a Riemann solver and the corresponding kinematic wave 
model can be used to study many other transportation network problems. 
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